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ABSTRACT 



Aims. We study the mean profiles of the multi-component pulsars PSRs B1839+09, B1916+14 and B2111+46. We estimate the 
emission height of the core components, and hence find the absolute emission altitudes corresponding to the conal components. 
Methods. By fitting Gaussians to the emission components, we determine the phase location of the component peaks. Our findings 
indicate that the emission beams of these pulsars have the nested core-cone structures. Based on the phase location of the component 
peaks, we estimate the aberration-retardation (A/R) phase shifts in the profiles. Due to the A/R phase shift, the peak of the core 
component in the intensity profile and the inflection point of the polarization angle swing are found to be symmetrically shifted in the 
opposite directions with respect to the meridional plane in such a way that the core shifts towards the leading side and the polarization 
angle inflection point towards the trailing side. 

Results. We have been able to locate the phase location of the meridional plane and to estimate the absolute emission altitude of 
both the core and the conal components relative to the neutron star center, using the exact expression for the A/R phase shift given by 
Gangadhara (2005). 

Key words. pulsar-PSRB 1839+09: B1916+14: B2111+46: Core emission height 



1. Introduction 

Pulsar radio emission is understood to be emitted by the rela- 
tivistic plasma accelerated along the dipolar magnetic field lines 
(e.g., Ruderman & Sutherland 1975). Among the various mod- 
els proposed for pulsar emission, the coherent curvature radia- 
tion has turned out to be an effective mechanism for explaining 
some of the important pulsar radiation properties. The common 
occurrence of an odd number of components in the mean pulsar 
profiles has lead to the nomenclature of a nested conal struc- 
ture for the pulsar emission beam (e.g., Rankin 1983a; Rankin 
1993). However, Lyne & Manchester (1988) suggested that the 
emission within the beam is patchy, i.e., the distribution of com- 
ponent locations within the beam is random rather than orga- 
nized in one or more hollow cones. Also studies by Mitra & 
Deshpande (1999) indicate that the structure of the pulsar emis- 
sion beam is more likely to be nested hollow cones. Gangadhara 
& Gupta (2001, hereafter GG01), and Gupta & Gangadhara 
(2003, hereafter GG03) showed that the prevalent picture of 
emission cones axially located around the central core compo- 
nent is a suitable model for explaining the core-cone structure of 
the pulsar emission beam. 

A long-standing question in pulsar astronomy has been the 
location of the radio emission region in the magnetosphere. In 
the literature, there are mainly two types of methods proposed 
for estimating the radio emission altitudes: (I) a purely geo- 
metric method, which assumes the pulse edge is emitted from 
the last open field lines (e.g., Cordes 1978; Gil & Kijak 1993; 
Kijak & Gil 2003), (2) a relativistic phase shift method, which 
assumes that the asymmetry in the conal components phase lo- 
cation relative to the core is due to the aberration-retardation 



phase shift (e.g., GG01, Gangadhara 2005, hereafter G05). Both 
methods have merits and demerits: the first method has an am- 
biguity in identifying the last open field lines, while the latter 
is restricted to the profiles in which the core-cone structure can 
be clearly identified. The emission heights of PSR B0329+54 
given in GG0 1 , six other pulsars in GG03 and the revised ones by 
Dyks, Rudak & Harding (2004, hereafter DRH04) are all relative 
to the emission height of the core, which is assumed to be zero. 
However, the core emission is believed to originate from lower 
altitudes than that of the conal components (e.g., Blaskiewicz 
et al. 1991; Rankin 1993). Hoensbroech & Xilouris (1997) esti- 
mated the emission heights at high frequency radio profiles for a 
set of pulsars. They suggested that the emission heights at high 
frequency can set an upper limit for the core emission height. 

By assuming a fixed emission altitude across the pulse, 
Blaskiewicz et al. (1991, hereafter BCW91) presented a rela- 
tivistic rotating vector model. The results of this purely geomet- 
ric method are found to be in rough agreement with those of 
BCW91. However, the relativistic phase shift method clearly in- 
dicates that the emission altitude across the pulse window is not 
constant (GG01; GG03; DRH04; Johnston & Weisberg 2006; 
Krzeszowski et al. 2009). 

By considering the relativistically beamed radio emission in 
the direction of the magnetic field line tangents, Gangadhara 
(2004, hereafter G04) solved the viewing geometry in an in- 
clined and slowly rotating dipole magnetic field. A more exact 
expression for the relativistic phase shift is given in (G05), which 
also includes the phase shift due to polar cap currents. In the 
present work, we analyze the mean profiles of PSRs B 1839+09 
and B1916+14 at 1418 MHz, and PSR B21 11+46 at 610 MHz 
and 1408 MHz, to estimate the absolute emission height of the 
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Fig. 1. Schematic diagram showing the A/R phase shift between the 
core peak (CP) and the polarization position angle inflection point 
(PPAIP). The panel (a) for the co-rotating frame, where the phases of 
PPAIP and CP coincide with that of the meridional plane (M), and (b) 
in laboratory frame, due to A/R effects, both the PPAIP and the CP are 
symmetrically shifted in the opposite directions with respect to M. 



pulse components. In Sect. 2, we give a method for estimating 
the absolute emission height of pulse components. 

2. Method for estimating the absolute emission 
heights 

The work of BCW91 generalized the rotating vector model 
(RVM) to include the relativistic effects due to rotation. 
According to their model, the centroid of the intensity profile 
advances to an earlier phase by ~ r/r^c, while the polarization 
position angle inflection point (PPAIP) is delayed to a later phase 
by ~ 3 r/r^c, where r is the radial distance from the center of 
neutron star and r^c is the light cylinder radius. After estimat- 
ing the width of the pulse at * 10% intensity level and by fitting 
the relativistic rotating vector model, they estimated the phase 
shift between the centroid of the profile and the PPAIP. But the 
retardation phase shift was ignored in BCW91, as they assumed 
a constant emission height across the pulse. In GG01, GG03, 
Johnston & Weisberg (2006) and Krzeszowski et al. (2009) it 
was shown though that the emission altitude is not constant 
across the pulse, and hence retardation has to be taken into ac- 
count while estimating the A/R phase shifts. Further, Dyks et al. 
(2004) showed that the centroid of the intensity profile advances 
by ~ 2 r/r^c while the PPAIP is delayed by ~ 2 r/r^c due to A/R 
effects with respect to the meridional plane. 

By solving the viewing geometry in the dipole magnetic 
field, Gangadhara (2005) showed that instead of a centroid of in- 
tensity profile, the phase shift of the central (core) peak (CP) rel- 
ative to the meridional plane must be considered for estimating 
the A/R phase shifts. In the co-rotating (non-rotating) case both 
core and PPAIP originate from the same phase (meridional plane 
M, see Fig. la). Whereas in the observer (laboratory) frame, the 
CP shifts to the earlier phase and the PPAIP to the later phase 
by the same magnitude (see Fig. lb). Hence to find the absolute 
emission height of the profile components including the central 




Fig. 2. Schematic diagram to show the probable distribution of emis- 
sion patterns across the pulsar beam. The beam cross sections as they 
appear in (a) the co-rotating frame and (b) the laboratory frame, where 
the cones are not coaxial with the central core because of A/R retar- 
dation effects. The thick horizontal line represents the direction of the 
tracing of the line-of-sight across the beam. The resultant intensity pro- 
file is shown in the adjoining box. The vertical line denotes the merid- 
ional plane, and the thick arrow represents the direction of the pulsar 
rotation. 



(core) component, we adapted the method of G05 to consider 
the CP instead of the centroid of pulse (BCW91) for estimating 
the A/R phase shift. It is logical to presume that the aforesaid 
r should be the same for the origin of the central (core) com- 
ponent and the PPAIP. Or stated otherwise, the phase difference 
A0£ P = (p' core - 0p PAIP ~ 0, in the co-rotating frame, where cf>' cole 
is the phase location of the core peak while PPAIP is the phase 
location of the PPAIP. 

As illustrated in the Fig. [2j the panel (a) depicts the cross 
section of the emission region in the co-rotating frame, while 
panel (b) shows for the same in the observer's frame. The thick 
arrow represents the direction of the rotation, and the thinner line 
the sweep of the line-of-sight. The shaded ring-like region rep- 
resents the nested conal emission regions and the central circle 
the core emission region. The sweep of the line-of-sight across 
the depicted region causes the core peak and the PPAIP to be 
separated by a roughly equal measure (~ 2r/rLc) in opposite di- 
rections from the meridional plane (M), as illustrated in Fig.Q] 
The resultant intensity profile, which characterizes a sum total 
of emissions after the line-of-sight sweeps across the emission 
region, is shown in the adjoining box on the right hand side. 

3. Application of the method 

The implementation of the aforesaid method to estimate 
the absolute emission height of the core and conal compo- 
nents is described below. We considered the mean profiles of 
PSRs B1839+09, B1916+14 and B21 11+46 for our study, as 
they exhibit a clearly identifiable core and smooth polarization- 
position-angle (PPA) swing. We obtained the data of PSRs 
B1839+09 and B1916+14 from Everette & Weisberg (2001), 
and those of PSR B2 1 1 1 +46 from EPN data base. 

3. 1 . Longitude of the core peak 

We fitted Gaussians to the pulse components to resolve the in- 
dividual components based on the method developed by Kramer 
et al. (1994), and the profiles are given in Figs. [3]— [6j Hence the 
peak-phase locations of the individual components are resolved. 
The broken line curves in panel (a) indicate the fitted Gaussians. 
The phase location of the core peak (<f>' C0le ) is marked with an 
arrow and is tabulated in Table Q] The PPA is plotted in panel 
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(b), and the vertical lines mark the fitted region of the curve. The 
panel (c) shows the zoomed-out region of the PPA within the re- 
gion of the fit. The arrow points to the PPAIP in both panels (b) 
and (c). 

The location of the central component (core) is expected to 
appear at M for an observer in the co-rotating frame as explained 
in Sect.|2]and illustrated in Fig. 1(a). But for an observer in the 
laboratory frame, the core emission will be advanced to an ear- 
lier phase by 6(/>' cole 2 r cole /rLc and the corresponding PPAIP 

delayed to a later phase by <*>0 PPAIP ~ 2 r core /rLc, where r core is 
the emission height of the core. Then the A/R phase shift of the 
core with respect to M is 6<f>' CQK = A<f>' /2, and the parameters re- 
lated to core emission are given in Table Q] The frequency v of 
each data set is given in col. 2, and the phase shifts 6(/>' cole and 
6<p' 



PPAIP 



in cols. 3 and 4. 



3.2. The relativistic RVM fit 

We fitted the relativistic RVM (BCW91) to the region of the 
PPA data, which corresponds to the core emission, the region 
over which the core polarization is significantly higher than that 
over the adjacent conal regions. We could see from the pro- 
files that the pulse phase range falling within the full-width-at- 
half-maximum (FWHM) of the central-fit-Gaussian can clearly 
bracket out the emissions that dominate the core. The core emis- 
sion is found to be dominant within 10% intensity levels of the 
central Guassian of the 1418 MHz data of PSR B1916+14, and 
hence we included that region for fitting the BCW91 curve for 
that profile. Additional reasons for restricting the fit region are 
discussed in detail in later sections. 



3.2.1. Longitude of polarization position angle inflection point 

We invoked the guess values for the emission height r and the 
phase 0p PAIP to fit the BCW91 curve. The polarization angle 
data, falling within the FWHM of the core, are fitted with a 
6th degree polynomial: iff (</>') = flo + ci\4>' + «20' 2 + • • •> where 
ao, a\, «2, • • • are the fit coefficients, and <p' is the pulse phase 
in degrees. We differentiated the fitted polynomial and found the 
maximum of \di[r/d(f>'\ that gives the guess value of PPAIP - In 
the next step, the PPA data were fitted with the relativistic RVM 
curve (BCW91) using the following expression 



<Abcw = iAo + tan 



sin a sin(fif) - 3(r/rix) sin^" 



sin ft + sin a cos ((1 - cos(£2f )) 



(1) 



where £ = a + ft. The inclination angle of the magnetic axis 
relative to rotation axis is a, and the impact angle of the line-of- 
sight relative to the magnetic axis is ft. The fitted PPA data are 
shown in the panel (b) of Figs. [3]-|6] 

The parameter </<o is inserted in the above Eq. ([TJ in order to 
offset the possible arbitrary and constant 'vertical shift' that the 
raw PPA data might have. This is due to the arbitrary value of 
the projection of the rotation axis in the sky plane. Since we are 
interested in finding r, which is expected to be relatively constant 
in the region of fit, we assumed it to be independent of 0' (= Q. t), 
and hence it is taken as a fitting parameter. The vertical shift in 
the raw PPA data, in the first step, was brought closer to zero by 
finding a trial value for the PPAIP on the vertical axis (i/^r) from 
the polynomial fit, and thereafter the data were shifted vertically 
so that i^pol — > <Apol _ iffj, where i/'pol is the observed polar- 
ization angle. These PPA data were fitted with the Eq. (Q]i keep- 
ing r and if/^ as the free parameters; hence allowing two degrees 
of freedom for the fit function, i.e., allowing the fit function to 
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Fig. 3. Intensity profile of PSR B1839+09 at 1418 MHz, fitted with 
the Gaussians to the sub-pulse components. In panel (a) the continuous 
line represents the observed mean profile while the broken line curves 
represent the fitted gaussians. The arrow points to the phase of the core 
peak. In panel (b) the corresponding polarization angle (x) is fitted with 
relativistic RVM curve (BCW91) curve. The arrow points to the phase 
of inflection point (PPAIP), and the vertical lines mark the region of 
phase over which the BCW91 curve is fitted. In panel (c) the zoomed 
out region demarcated by the vertical lines in panel (b), is plotted and 
the arrow points to the PPAIP. 



'adjust' in both the vertical (through the free parameter and 
horizontal direction (through free parameter r) for a good fit. 
The fit procedure was repeated consecutively a few times, with 
the values of r and ifto from the preceding fit as guess values for 
the final convergent and stable values. Thus PPAIP was found 
from the fit, and the corresponding values are shown in Table [T] 
The geometric angles a and ft were not invoked as fit parame- 
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ters in Eq. ([TJ because we used their predetermined (published) 
values in the BCW91 fit function. We used the values of a and 
P given by Everette & Weisberg (2001) for PSRs B1839+09 and 
B1916+14, and for PSR B21 1 1+46 by Mitra & Li (2004, here- 
after ML04). 



3.3. The longitude of the meridional plane 

The CP and PPAIP of the polarization angle appear at M in the 
co-rotating frame, as indicated by Fig. QI a )- But they are sym- 
metrically shifted in opposite directions from M due to the A/R 
phase shift for an observer in the laboratory frame, as indicated 
by Fig. |TJb). Note that the phase location of M is invariant with 
respect to the rotation effects, or in other words, both the CP and 
the PPAIP come closer to M at smaller r, and move away from 
it at larger r. If (f>' COIf , and <p' PPA1P are the estimated phases of the 
CP and the PPAIP, then their phase difference (A</>') is given by 
A</>' = PPAI p-0core' an d the meridional plane M is situated at the 
mid point between CP and PPAIP. Hence the phase of M is given 
by <p' M = <p' core + (A</>'/2). In Figs.|3]-|6]the phases are shifted by 
(f>' M , so that M appears at the zero phase and an arrow points to 
the phase location PPAIP in the polarization angle panels (b). 



3.4. The A/R phase shift of the core 

The location of the central component (core) should appear at M 
for an observer in the co-rotating frame as explained in Sect|2] 
But for an observer in the laboratory frame, the core emission 
will be advanced to an earlier phase by S<p' colf , ~ -2r core /rLc 
and the corresponding PPAIP be delayed to a later phase by 
<S0 ppAI p ~ 2r core /rLc, where r core is the emission height of the 
core. Then the A/R phase shift of the core with respect to M is 
<%, re = A</>'/2. 



3.5. Phase locations of the core and cone component peaks 

We found the peak locations of the individual components by 
fitting Gaussians to the mean profiles of the three pulsars. As 
mentioned in Sect. 13.31 the data were shifted by <p' M so that 
M appears at zero phase. Because the A/R effects are absent 
in the co-rotating frame, the conal components are expected to 
be symmetrically located on either sides of meridional plane as 
indicated by Fig. [3a). But in the laboratory frame, the cones 
are advanced to an earlier phase due to the A/R effects and 
are hence asymmetric with respect to the core location as indi- 
cated by Fig. EJb). The meridional plane M is taken to be at the 
zero phase, and the measured phases are, therefore, the absolute 
phases with respect to M. Accordingly the estimated emission 
heights are the absolute emission altitudes with respect to the 
center of the neutron star. 

Let <p' h and (f>' T be the peak locations of the conal components 
on the leading and trailing sides of a pulse profile, respectively. 
Then, using the following equations, we estimate the A/R phase 
shift (£</>') of the cone center with respect to M, and the phase 
location (</>') of the component peaks in the absence of the A/R 
phase shift, i.e., in the co-rotating frame (see Eq. (1 1) in GG01): 



Sep' = i(_4>' T + cf > ' I ), 0' = I(^-^). (2) 
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Fig. 4. Intensity profile of PSR B1916+14 at 1418 MHz. See the cap- 
tion of Fig.f3]for details. 

4. The absolute emission heights 

4.1. Emission height of the core 

The core emission height was computed by using the 6<p' C0[e in the 
expression for the A/R phase shift given by G05 (see Eq. (45)). 
The parameters related to the core emission are given in TableQ] 
In Col. 6 the emission heights are given as a percentage of the 
light cylinder radius tlc- It shows that the core emission in the 
radio band occurs over a range of altitude spanning from 0.2 to 
1 per cent of the light cylinder radius. The radio frequency v of 
each data set is given in Col. 2, and the phase shifts 6<p' core and 
^ppaip m Cols. 3 and 4, respectively. The values of x 2 and the 
standard residuals obtained are given in Cols. 7 and 8, respec- 
tively. The foot location of magnetic field lines on the polar cap 
relative the magnetic axis are given in Col. 9. 
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Table 1. Core emission geometry parameters of PSRs B1839+09, B1916+14 and B21 1 1+46 



5 



Pulsar 


V 


A/A' 


A/A' 

^PPAIP 


'"core 


'"core 


X 






PSRB 


(MHz) 


(°) 


C) 


(km)" 


(% 1c) 




(%? 


foot value 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1839+09 


1418 


-0.32±0.07 


0.32±0.01 


50+5 


0.28+0.03 


1.02 


100 


0.51+0.03 


1916+14 


1418 


-0.23±0.02 


0.23+0.01 


95+4 


0.20+0.01 


28.25 


94 


0.26+0.01 


2111+46 


610 


-1.33±0.02 


1.33+0.20 


503+38 


1.04+0.08 


5.7 


100 


0.13+0.00 




1408 


-0.23±0.08 


0.23+0.29 


83+54 


0.17+0.11 


2.05 


95 


0.51+0.17 



"Emission heights computed using the exact formula (G05). 

*The percentage of the standardized residuals (SR) with value of -2 and +2. 

Table 2. Conal emission geometry parameters of PSRs B1839+09, B1916+14 and B21 1 1+46 



Pulsar 


V 


Cone 




<t>' t 


6<p' 


r 


r 


r 


■sMof 


PSRB 


(MHz) 


(No.) 


o 


O 


C) 


n 


(km)" 


(% r LC ) 


foot value 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1839+09 


1418 


1 


-3.75+0.15 


2.96+0.18 


-0.39+0.12 


4.06+2.30 


63+19 


0.34+0.10 


0.80+0.10 


1916+14 


1418 


1 


-3.49+0.05 


2.04+0.07 


-0.72+0.04 


2.63+1.00 


299+18 


0.63+0.04 


0.39+0.01 


2111+46 


610 


1 


-14.78+0.45 


11.79+0.31 


-1.49+0.27 


2.85+1.10 


581+107 


1.20+0.22 


0.30+0.02 






2 


-33.83+0.21 


24.95+0.22 


-4.44+0.15 


5.87+1.10 


1891+64 


3.91+0.13 


0.35+0.00 




1408 


1 


-14.70+0.40 


11.78+0.34 


-1.46+0.27 


3.30+1.80 


533+97 


1.10+0.20 


0.37+0.03 






2 


-28.47+0.31 


22.48+0.34 


-2.99+0.23 


5.58+1.80 


1151+87 


2.38+0.18 


0.42+0.01 



"Emission heights computed using the exact formula (G05). 



4.2. Emission height of the cones 



We found the emission height of the cones based on the pro- 
cedure that is described in Sect. 13.51 In Col. 3 of Table [2] we 
have given the cone numbers and the peak locations of the conal 
components on the leading and trailing sides in Cols. 4 and 5, 
respectively. The conal components are believed to arise from 
the nested conal emissions (Rankin 1983a, 1983b, 1990, 1993), 
which along with the central core emission make up the pulsar 
emission beam. In Col. 6 of Table [2] we have given the values 
of dip' . In general it increases in magnitude from the innermost 
cone to the outer cone. The half-opening angle T (see Eq. (7) 
in G04) of the emission beam is given in Col. 7. Using the ex- 
act formalism for the A/R phase shift (see Eq. (45) in G05), we 
computed the emission heights and show them in Col. 8. Their 
percentage values in tlc are given in Col. 9. In Col. 10, we give 
the normalized co-latitude s/s\ f of the foot field lines on the 
polar cap, which are associated with the component emissions. 
Due to the relativistic beaming and restrictions owing to geome- 
try, we find that the observer tends to receive the emissions from 
open field lines, which are located in the foot locations ranging 
from (approximately) 0.1 (Table 1) to 0.8 (Table 2) on the polar 
cap. 



4.3. PSRB1 839+09 

By fitting three Gaussians to the mean intensity profile of PSR 
B 1839+09 we have identified three sub-pulse components: a 
central component and a pair of outer components flanking the 
central one. By invoking the picture of a nested cone structure 
we infer that the outer pair of components corresponds to a conal 
emission. The inner component might be due to emissions close 
to the magnetic axis or because of the 'grazing-cut' of the line- 
of-sight with an inner ring of emission. The point of emission 
for the central component should fall in the meridional plane 
M in the co-rotating frame of the pulsar in either of the cases. 
We clearly identified the phase locations for the PPAIP and the 
CP relative to the meridional plane M for PSR B 1839+09. The 
absolute emission heights estimated for the core and conal com- 
ponent are given in Tables [TJ and [2] respectively. The emission 
heights of the central (core) component and cone are found to be 
~ 50 km and ~ 60 km, respectively. 

4.4. PSRB1916+14 

By fitting Gaussians to the intensity profile, we could identify 
three sub-pulse components: a central component and a pair of 
outer components flanking the central one. We clearly identified 
the phase locations of the PPAIP and the central peak, and found 
the phase of meridional plane M. The absolute emission heights, 
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Fig. 5. Intensity profile of PSRB21 11+46 at 610 MHz. See the caption Fig. 6. Intensity profile of PSR B21 11+46 at 1408 MHz. See the cap- 



of Fig.[3]for details. 



tion of Fig.[3]for details. 



estimated for the central (core) and conal components, are given 
in Tables [T]& |2] The emission heights of the central component 
and the cone are estimated to be ~ 100 km and ~ 300 km, re- 
spectively. 

4.5. PSRB2111+46 

It is a well-studied pulsar. By fitting Gaussians to its intensity 
profiles at frequencies 610 MHz and 1408 MHz, we could iden- 
tify five sub-components: a central (core) component, a pair of 
inner components and a pair of outer components. One can guess 
at the presence of two inner components hidden between the core 
and the outer components even by visual inspection. These two 
hidden components were detected in the previous work of GG03 
in the 333 MHz single pulse data. The absolute emission heights 



estimated for the core and the cones are given in Tables Q] and [2] 
The core emission height is found to be ~ 500 km at 610 MHz 
and ~ 80 km at 1408 MHz. It has been reported in ML04 that 
the widths of the core tend to show a significant frequency evo- 
lution in their chosen set of six pulsars, and hence they argued 
that the core emission does not come from the stellar surface. 
However, we need to consider more high quality data sets at 
different frequencies to see the frequency evolution of the core 
emission height, and whether it follows any radius-to-frequency 
mapping. 



5. Results and discussions 

Based on the A/R method, we estimated the absolute emis- 
sion height of the core as well as cones in three pulsars: PSRs 
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B1839+09, B1916+14 and B21 11+46. Though this method is 
based on the existing standard models in literature, the combi- 
nation of the A/R phase shift and the delay-radius relation of 
BCW91 for estimating the core height is novel. 

The geometrical method, involving a comparison of the mea- 
sured pulse widths with geometrical predictions from dipolar 
models, is believed to yield absolute emission heights. However, 
the estimation of emission height, using the geometrical method, 
is based on the assumption that the pulse edges originate from 
the last open field lines of the polar cap. In general, the edge of 
the on-pulse region may not originate from the last open field 
line, and hence the assigning of the edges of the intensity profile 
to the last open field lines can be misleading. For example, the 
range of magnetic foot-colatitude for field-lines that are associ- 
ated with components in PSR B21 1 1+46 are in the range from 
S/Siaf ~ 0.13 to 0.5, whereas the last open field line is at 1. This 
means that the boundary of the active region of emission can lie 
anywhere from ss 0.5 to 1 . 

According to DRH04, the A/R phase shift advances the 
centroid of the intensity profile to an earlier phase by 5<p' c = 
2nof/'"LC, while the PPAIP is delayed to a later phase by 
(50p PAIp ~ 2r core /rLc, where r\ f is the emission height from 
the last open field-line and r core is the emission height of the 
core. Then A<p' = 2(n f + r core )/rLc, and the emission height 
r - tlc A<j>'/4 = (nof + >"core)/2, gives only an average of the 
emission height for the core and the pulse edge, which is far from 
the true value. This emission height cannot represent any spe- 
cific pulse sub-component of the profile, and can be misleading 
in cases where r\ Q f and r core are significantly different. Further 
more, this will introduce large systematic errors in the emission 
heights estimated from geometrical methods, due to the afore- 
said assumption of identifying the last open field lines with the 
pulse edges. 

Rankin (1983a) has argued that the pulsar emission cones are 
quasi-axial, i.e., the conal components are not exactly axially 
located with respect to the magnetic axis. Mitra & Deshpande 
(1999) have suggested that the pulsar emission beams are nearly 
circular in the aligned configuration (a ~ 0°) and change to ellip- 
tical in the orthogonal configuration (a ~ 90°). The majority of 
the pulsar observations indicate that the beam geometry is likely 
to be nested cones, distributed in a nearly non-coaxial fashion 
about the magnetic axis. A likely case is that the cones, which 
are coaxial in the co-rotating frame, will appear non-coaxial in 
the laboratory frame because of the A/R phase shifts (GG01; 
GG03). 

In the works GG01 and GG03, the emission height of the 
core was neglected by assuming that it is considerably smaller 
than that of the components. However, we find that the emission 
height of the core is quite significant and cannot be neglected in 
comparison to the emission height of the components. We iden- 
tify the meridional plane M as being located at the mid point 
between the centroid of the intensity profile and the PPAIP, ow- 
ing to the A/R effects. By recognizing this, we were able to esti- 
mate the absolute emission heights of both the core and the conal 
components. 

As mentioned before, we restricted the region of the fit of 
the BCW91 (relativistic RVM) curve to the section of the PPA 
data falling within the FWHM of the core component for esti- 
mating the core emission height, and the justification for doing 
so is given now. The expression for the BCW91 was derived 
by assuming that the emission altitude across the active region 
of the pulse profile is a constant. Thus in a BCW91 fitting, a 
single r value was taken to characterize the emission height of 
the full region of the PPA data. But later observational results 



(e.g., GG01, GG03) established that the emission altitude corre- 
sponding to the subpulse components in multi-component pro- 
files spans over a large range of emission heights. This elicits 
the fact that in multi-component profiles the r, found by fitting 
the BCW91 curve to the full PPA region of the active profile, 
might give an emission altitude that can be significantly differ- 
ent from those obtained from the A/R method for the subpulse 
components. 

The best-fit value of r in the BCW91 model, which is the 
weighted average of r; that characterize the emission height at 
each point of pulse phase, is given by Eq. ( IC.8b in Appendix C 



Hence a single value of r found from the BCW91 fit cannot be 
closer to the true emission height corresponding to the core or 
cone peak if r\ varies significantly within the region of the fit. For 
example, one can compare the emission altitudes in our Tables Q] 
and|2]with those given in Table 3 of ML04 for PSR B2111+46. 
It can be surmised that a single value of r cannot characterize 
the emission altitude across the entire active region of a multi- 
component profile. 

We can think of two viable alternatives in this scenario: ei- 
ther (1) adapt or modify the BCW91 formulation for a variable 
emission altitude r (Dyks 2008) or (2) fit the BCW91 curve for 
regions of the PPA profile having a relatively constant value of r. 
We prefer the latter alternative to evade the modification of the 
theory behind the BCW91 model. Here we note that a and /3 are 
not invoked as fit parameters; instead we used their published 
values in Eq. (IC.ll l. This is expected to further reduce the am- 
biguity of the fit results and to aid in counteracting an obvious 
disadvantage in this 'restricted' fit method, i.e., having a reduced 
number of fitted PPA data points than a fit for the 'full range' of 
PPA data. It is remarkable that some of the fit statistics (e.g., 
reduced^ 2 ) given in Table Q] reveal that the present method of 
fitting is comparable (in a few cases even better) to the existing 
ones in the literature (e.g., compare the^ 2 value given in TableQ] 
with Table 2 of ML04 for PSR B21 11+46). We estimated the 
standardized residuals (SR) and found the percentage of SR that 
falls within -2 and +2 as given in Col. (8) of Table Q] As it is 
known, a good fit is expected to have a threshold 95 % of the 
standardized residuals to fall within -2 and +2. Several previous 
works found that a and /? are highly covariant in PPA fits (e.g., 
Everette & Weisberg 2001). But this covariance of a and /3 with 
the r parameter was not mentioned by any of them. The fit statis- 
tics do not reveal any significant covariance of a and /3 with r. 
This gives us a further clue for finding the r parameter without 
invoking a concurrent fit for a and /3 (see Appendix A for the 
fitting procedure). 

Owing to the extreme difficulties encountered in determining 
a and (3 through RVM fitting, a larger range of PPA data have al- 
ways been preferred for a better fit (e.g., Everette & Weisberg 
2001). The justification for doing so is that a and /3 must remain 
constant throughout the entire PPA profile. But in the present 
scenario, as described earlier, the selection of a large range of 
PPA data for fitting does not always translate into a better esti- 
mation of r because of the variation of the emission height with 
pulse phase. So, owing to all of the above said reasons we restrict 
the fit of the BCW91 curve to the PPA profile, falling around the 
core component, which is expected to yield an emission altitude 
characterizing the core height. 

A section of inner cones often lapses over the core as is seen 
in the Gaussian fits (panel (a) of Figs. [3}j6]l of the total inten- 
sity profiles. The inner cones may contribute to the core polar- 
ization near the edges of pulse phase of the FWHM region that 
we bracketed. Hence the PPA corresponding to the bracketed re- 
gion will be 'contaminated ' by the adjacent conals, and this has 



8 



Thomas and Gangadhara: Absolute Emission Altitude of Pulsars 



to be accounted for. We estimated and accounted for the error 
induced because of this effect in the estimation of the core emis- 



sion heights (see Appendix B and Appendix C I. 

The possibility that the A/R phase shift may be reduced 
by the rotational distortion of the magnetic field line due to a 
sweep-back of the vacuum dipole magnetic field lines has to 
be considered. The sweep-back of dipole magnetic field lines 
was first treated in detail by Shitov (1983). Further, Dyks & 
Harding (2004) investigated the rotational distortion of pulsar 
magnetic field by making the approximation of a vacuum mag- 
netosphere. For <p' = 30°, /3 = -1.6° and a = 14° we com- 
puted the phase shift S(p' m{sh due to the magnetic field sweep-back 
(Dyks & Harding 2004^0 see Eq. (49) in G05). It is found to 
be < 0.0001 rad for r/r^c < 0.06, which is much smaller than 
the aberration, retardation and polar cap current phase shifts in 
PSR B21 11+46. Hence we neglect the magnetic field sweep- 
back effect. 

The field-aligned polar-cap current does not introduce any 
significant phase shift into the phase of the PPAIP. But it intro- 
duces a positive offset into the PPA, though it roughly cancels 
due to the negative offset by aberration (Hibschman & Arons 
2001). The phase shift of pulse components due to the polar cap 
current was estimated recently by G05, and found to be quite 
small compared to the A/R phase shift. 
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6. Summary 

We analyzed the mean profiles of PSRs B1839+09 and 
B1916+14 at 1418 MHz, and those of B21 11+46 at 610 MHz 
and 1408 MHz. The phase of the peak of central component 
(core) and that of the polarization position angle inflection point 
are symmetrically shifted in the opposite directions with respect 
to the meridional plane due to A/R effects. By recognizing this, 
we were able to locate the phase of the meridional plane and 
to estimate the absolute emission altitudes of the core and the 
conal components relative to the center of the neutron star. We 
used the exact expression for the phase shift given recently by 
G05. In all the cases we found that the core emission occurs at a 
relatively lower altitudes than the conal emissions. It is also in- 
teresting to note that the core emission at different frequencies in 
PSR B2 1 1 1 +46 falls in a range of altitude of 80 km at 1408 MHz 
to about 500 km at 610 MHz. It is clear that the low frequency 
emission comes from a higher height than that at high frequency. 
However, to confirm whether the core emission heights also obey 
any radius-to-frequency mapping demands the recursive analy- 
sis with high quality multi-frequency data. We plan to employ 
the methods described in this paper for the study of a few other 
pulsars with high quality data. 
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Appendix A: Estimation of 1 cry 

We apply the standard methods of statistics for fitting the ob- 
served data with the model. For fitting the PPA data with the 
BCW91 curve we define the reduced^- 2 as 



X 



-s 



lApOlM) - ^BCW(0,O 



(A.l) 



where *Ppol is the observed polarization angle and ^bcw is the 
model (BCW91) value at the discrete rotation phase <p' x , and N is 
the number of data points. Minimization of x 1 gives the best fit. 
We use Mathematica version 7 for the fitting and other calcula- 
tions that ensue therewith. 

For profile regions with a very high value of L/crj, i.e., for 
L/cri » 10, the cr^j is taken as 



o> = 



2L 



(A.2) 



where ctq and <t\j are the RMS values of the off-pulse Q and U 

stokes parameters, respectively. The parameter L = ^Q 2 + U 2 
represents the linear polarization and o-j is the standard devia- 
tion of total intensity in the off-pulse region. Though the mea- 
surements of U and Q around their true values are normally dis- 
tributed, the PPA measurements are not similarly distributed at 
the intermediate values of L/crj. Hence the form of cr^ as given 
above is not appropriate in this regime. The appropriate distribu- 
tion function that characterizes the ifr distribution is discussed in 
detail by Naghizadeh-Khouei & Clarke (1993). The application 
of their scheme is discussed by Everette & Weisberg (2001); we 
follow their method for finding the <r^ for the PPA data points. 
The probability distribution function of the measured position 
angle if/ around a true value i/f true is given by Naghizadeh-Khouei 
& Clarke (1993) as 



GOA; <Ao, Po) = -4= (-7= + W^O + erf(T] ))\e- 



(A3) 



where 770 = (P / V2) cos[2((/r - i/w)], Po = LtmeM, erf is 
the error function and Lt me is the unbiased linear polarization 
found from the measured linear polarization. Hence the 1 cr^ con- 
fidence level in the PPA is found out by adjusting the limits of 
integration and fixing if/ tme = 0, 



f 

J-i. 



G((A;A))# = 62.86%. 



(A.4) 



We set up a table of the integration values of the integral against 
a series of discrete Po values and the lcr^ levels at integral value 
of 0.6286 are found by interpolation. 

Appendix B: Polarization angle *F: Contribution 
from the adjacent component to the core 

The polarization angle ¥ is defined as 
2 \Q 



Defining 

U = U Q + Ui , 

Q = Go + Qi , 
l = Vt/ 2 + Q 2 , 

L0 = yjlJl + Ql , 



(B.l) 

(B.2) 
(B.3) 
(B.4) 

(B.5) 



where U\ and Q\ are taken as sufficiently small, ¥ can be a series 
expanded up to the first order in U\ and Q\ as 



To = -tan 

2 



2o 



2 f/ ( 2 + e ( 2 , 



The expression for A^V may be approximated as 

a I QU1-UQ1 
2 L 2 



(B.6) 
(B.7) 

(B.8) 



(B.9) 



provided Ui and Qi are sufficiently small, so that Lq L. 

We use the above expressions for estimating the contribution 
to T from the adjacent conal components. The suffix "l" indi- 
cates the U and 2 contribution solely from the inner cone, while 
the suffix'0' indicates the pure core contribution for the same. 
We employ approximations to find the value of U\ and Q\ . 

Since the total intensity of an inner conal component can be 
fitted with a Gaussian, we make an assumption that the U\(<p') 
and Qi((f>') also follow a Gaussian shape within the adjacent 
conal component. Thus 



C/i(0') = f/l(max)exp 



and 



2l(0') = 2l(max) exp 



2\ 



(B.10) 



(#nax - 0') 



2\ 



(J A. 



(B.ll) 



where cr^ is the width and <p' m&x is the peak phase of the inner 
conal component (same as the corresponding Gaussian for the 
total intensity), while f/i(max) and Qi(max) are the peak values 
of U\ and Q\, respectively. The total intensity of the inner cone 
at phase <f>' i v where the Gaussians corresponding to the core and 
outer cone intersect, will have the minimum. Hence the values of 
Ui(4>' int ) and Qi((f>' in[ ) may represent approximately the true val- 
ues of U 1 and Q\ for the inner-cone at i 
and ( 1B.1U we can write 



£/l(max) = Ui((f>[ at ) exp 
Gl(max) = Ql(4>[J exp 



(0max - 4>[J 



(0max - 0[ nt ) 



2^ 



. By using Eqs. ( IB. 10b 

(B.12) 
(B.13) 



Using the Eqs. dB~T2l and (lB~T3l . we can find U x ((p') and Q\((p'), 
and hence the A v F(0') is estimated within the bracketed region of 
PPA. 

Note: The PSR B 1839+09 has practically no contribution of 
polarization from the adjacent conals to the bracketed core re- 
gion. Hence this analysis is not performed for it. The £/i(max) 
and 2i( max ) f° r PSR B1916+14 are the values at the peak 
phases of the adjacent cones and are directly found from the 
profile data. So, the Eqs. ( IB. 12b and ( IB. 13b are not applicable 
for it. Hence the above-said analysis is done for the profiles of 
PSR B21 11+46 only. 
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Appendix C: Determination of the r parameter: 
Correcting for adjacent conal contribution 

Here, we explain the scheme of the determination of r and 
^ppaip usm 8 the BCW91 method with predetermined values of 
a and ft. Consider the expression for the polarization angle given 
inBCW91: 



0BCw(0i) = tan 



Here define x 2 as 



sin or sin0; - 3(r/rLc) sin£ 



sin p + sin a cos £( 1 - cos 0;) 



X 



z 



tan[ PO l(0O] - tan[ B cw(0P] 



0"tan(*(0P 



(CI) 



(C2) 



where cr tan ^(0p = tan[cr^(0p], cr^(0p is the lcr error deter- 
mined by using the Naghizadeh-Khouei & Clarke ( 1 993) method 
( | Appendix A) . The definition of x 1 as given in Eq. ( 1C.2I > is 
slightly different from that used in the fitting, and this is meant 
exclusively for an explanatory purpose. Using the^- 2 of this form 
makes it easy to derive the following expressions for the r param- 
eter from the BCW91 model. The form of x 1 that is used in the 



actual fitting is given in Appendix A (Eq. IA.lt . and the method 
of correction for the adjacent conal contribution is briefly given 
at the end of this section. 

By minimizing the^ 2 with respect to the parameter r, we can 
write 



i 2 N 

d_X_ = 2 y 
dr t-L 

I 

d 



dr 



tan( iApolOP) - tan( B cw(0P) 
tan(0 B cw(0p) = 0. 



By substituting for tan(2 B cw(0OX we obtain 



3 

where 
Fi{<P[) -- 
and 

F(4>[) = 



r LC y FWmM i" sinf F 2 ^') 



tan t// 



(#) 



tan ifr 



(C3) 



sin /3 + sin a cos £( 1 - cos 00 



1 



sin p + sin a cos 1 - cos 00 



tan(0 POL (0 1 ')) (C.4) 



(C5) 



If r\ is a value so that 0pol(0i) = 0bcw(0O, (n may charac- 
terize the true emission height at the phase 0j), then by taking 
the analogy from Eq. dC.lt we can equate 



^i(0O 



3n 



sin f 



r L c sinyS + sinacos^(l - cos0p 
Then by substituting Eq. ( IC.6t into Eq. ( lC.3t . we get 



z< 



^ 2 (0P 



Z 



^ 2 (0O 



tan if/^ 



(C6) 



(C7) 



Within the bracketed region the small angle approximation for 
0' can be used so that we can re-write expression dC.7t as 
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1 



tan (// 
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(C8) 



or, by using S<f>' PPAlP 
write 



2r/r hC and 



J \ (shift) 



"i (shift) 2 , 



z 



tan ifc 



(0,0 



2rj/rLc we can 



(C.9) 



These expressions ( lC.8t and dC.9t show that the BCW fit of the 
PPA data is a weighted average of the emission heights r; of 
each data point. If some of the data points are corrupted by the 
adjacent conals, then the true linear polarization due to the core 
alone should be less than the observed. Hence it implies that 
the true PPA value corresponding to the pure core contribution 
should be slightly different from the observed PPA, which has 
a small mixture of contribution from adjacent conal component. 
An approximate method to estimate this error in the PPA (AT) 
due to the 'contamination' of the adjacent core is explained in a 
previous section. Hence the new weight factors are to be defined 
as 

CTTOTAL(0O = A / (r m„/^) +tan2 [ A>I, co„e(0O] , 

which takes into account the 'contamination' of the core due to 
the adjacent cones. Thus we find 



pure 



and 
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TOTAL 
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, /pure 
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"shift 2 
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TOTAL 



TOTAL 



Z' 



TOTAL v 



(CIO) 



(C.ll) 



using the expressions dC.8t and (IC.9t . In principle, r pme and 
^pfaip characterize the values of emission height and the cor- 
responding phase shift of the PPAIP found for the core region 
after correcting for the conal contamination. The error induced 
in locating the PPAIP due to the adjacent conal contamination 
can be given as A<50p PAIp = \S<ppp AlP - <50 PPAIP I- Hence the im- 
proved value of the phase location of PPAIP can be expressed 
as 



^ppaip ± A /Error((50j 



PPAIP 



fit 



Ac50' 2 pAIp 



(C.12) 



tan i/A 



where Error(£0 ppAIp ) nt is the lcr confidence interval for the 
<50 ppAIp obtained from the BCW91 fit without invoking the cor- 
rection for conal contamination. 

In the actual fitting procedure we used the form of y 2 as 
given in Appendix A We estimated |A<S0 PPAIP | for the profiles, 
and the values of <J0 PPAIP (Table 1) were attributed with an er- 
ror factor as given by Eq. dC.12t . A set of weights are found in 

the form of ct TO tal = ^cr 2 + A*Fj? one . The BCW91 model is fit- 
ted within the bracketed region after ( 1 ) weighting the data with 
weights cr^j and (2) again by weighting the data with weights 
0" total- The case (1) will yield the phase shift for the PPAIP, 
while the case (2) should yield the phase shift for the PPAIP with 
reduced weights to the PPA points where the conal contributions 
are present. The difference in the phase shifts found by case (1) 
and (2) should characterize the extra increment (decrement) in 
core emission height due to the adjacent conal contributions. The 
square of this phase shift difference is added with the squared er- 
ror of the PPAIP phase shift, and the square root of this sum will 
give the improved error factor for the core emission height. This 
improved error factor will take into account the error induced in 
the estimation of the PPAIP in the bracketed region due to the 
adjacent conal contribution. 



